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0\ • Abstract 

We express the complex potential £ and the metrical fields lo and 7 of all 

r^^ , stationary axisymmetric vacuum spacetimes that result from the application 

of two successive quadruple-Neugebauer (or two double-Harrison) transfor- 
(N 
> 

C^ , which are in turn easily expressed in terms of multipole moments. Moreover, 

^^ ! we suggest how, in future papers, we shall apply our approach to do the same 

■^ 

7' 

U 

bX). 



X 



mations to Minkowski space in terms of data specified on the symmetry axis. 



thing for those vacuum solutions that arise from the application of more than 

two successive transformations, and for those electrovac solutions that have 

axis data similar to that of the vacuum solutions of the Neugebauer family. 
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I. INTRODUCTION 

In recent years considerable interest has been displayed in stationary axisymmetric so- 
lutions of the Einstein-Maxwell equations that are characterized by axis data of the form 

with 

n 

U = J2 Uaz""-" , (1.2a) 

a=0 

n 

^ = E ^-^""" ' (l-2b) 

a=l 
n 

W = Y. Waz""-" , (1.2c) 

a=X 

Where E and $ are the complex potentials of Ernst,0 z is the Weyl canonical coordinate, 
and the coefficients in the polynomials are complex constants. In a long series of papers, 
V. Manko and his associates^ have been evaluating the complex potentials and metric fields 
for particular assignments of the axis data. For each such assignment, they solve anew 
SibgatuUin's integral equation formulationB of a Riemann-Hilbert problem. The question 
naturally arises "Could not all these solutions be obtained at once, rather than in the piece- 
meal manner employed by Manko et al.?" Our objective, which will be partially achieved in 
the present paper, is to express the complex potentials and metric fields of all such solutions 
in terms of arbitrarily prescribed axis data of the form indicated above. 

To address this question in a systematic way, we shall divide the solution of the problem 
into three parts: 

1. The general solution of the n = 2 vacuum problem {V = 0). 

2. The general solution of the n = 2 electrovac problem. 

3. The general solution of the problem for all values of n. 



This procedure will enable us to illustrate the basic ideas within the simplest context (1), 
after which we shall assemble heavier artillery to cope with problems (2) and (3). We already 
know that it will not be possible to solve problem (3) completely without resorting to some 
numerical work, but the situation is not quite as grim as one might suppose. As we solve 
problems (1) and (2), we shall point out how we intend to extend the procedures that we 
have used there to the case n > 2. 

The vacuum solutions, which are the subject of the present paper, all belong to the 
Neugebauer family,Q i.e., they can be generated from Minkowski space through the action of 
n successive quadruple- Neugebauer Backlund transformations. Alternatively, these solutions 
can be constructed using n double-Harrison Backlund transformations,^ or that Kinnersley- 
Chitre transformationQ that corresponds to the latter Backlund transformation. While the 
solutions have been known for a long time, this is the first attempt of which we are aware 
to express everything directly in terms of arbitrarily specified axis data in the manner of 
Manko et al. 

The complex potential S of the solution that results from applying a succession of n 
quadruple-Neugebauer Backlund transformations to Minkowski space is given by 



£ 



where U is the 2n x 2n determinant 



U 



in which occur the 2x2 submatrices 
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'z-KaY + P' 
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The 2n x 2n determinant W is constructed from U by replacing the (2n — l)-st row of the 
latter determinant by K^ ■ ■ ■ K^^- It is left for the reader to verify that ^ := U/W is a. 
solution of the complex potential field equationQ 

(a* - i)v'e = 2rve ■ ve , (1.7) 

if the constants Ka {a = 1, . . . ,n) are either real or occur in complex conjugate pairs and 
the constants Xa (a = 1, . . . , n) satisfy 

XaX; = 1 when Ka = K* . (1.8) 

The Kinnersley-Chitre transformations that correspond to various combinations of real i^'s 
and complex conjugate pairs of i^'s can be effectively unified into a single complexified K-C 
transformation in which parameters that are real in the case of an ordinary K-C transfor- 
mation are replaced by complex valued parameters. This approach was first explored by 
Hauser, who showed that one can solve the Hauser- Ernst homogeneous Hilbert problem just 
as easily for members of the group SL(2, C) as for members of SU{1^ 1).S Rather than think 
of every possible partition of if' s into real ones and complex conjugate pairs as comprising 
a different family of solutions, it is more natural and a lot more convenient to consider 
these as real cross sections of one family of complexified spacetimesp which we shall dub the 
Neugebauer family, honoring the person who pioneered the systematic study of members of 
this family and to whom the determinantal expressions for the S potential are due. 

The field U is homogeneous of degree n and the field W is homogeneous of degree n — 1 
in the r's. From this it follows, in particular, that the spacetime is asymptotically flat, with 
a possible NUT parameter. On the symmetry axis, where p = 0, one has r^ = K^ — z, so U 
and W are, respectively, polynomials of degree n and n — 1: 



U = Y.UaZ^-\ (1.9a) 

0=0 
n 

W = Y.Waz''-'' . (1.9b) 

a=l 

In particular, Uq is the determinant that is constructed from U by the substitution 



-1 (a = 1, 2, . . . ,n) . 



:i.io) 



When Uq ^ one can, if one wishes, readjust the common factor in U and W so that 
Uq = 1. We shall refer to the resulting constants Ua, Wa (a = 1, 2, . . . , ra) as the axis data, the 
specification of which uniquely determines the stationary axisymmetric vacuum spacetime. 
When Uq = 1, the real and imaginary parts of the constants Ua,Wa, {a = 1,2, ... ,n) 
are closely connected with the multipole momentsu'E^ that describe this asymptotically flat 
spacetime. If one translates along the z-axis so that Re Ui = 0, then iUi is a rotation 
parameter, and Wi a (complex) mass parameter, and so on. The imaginary part of the 
latter parameter is associated with the so-called NUT parameter of the spacetime. 

In general U and W can be expressed (up to a common constant factor) in the respective 
forms 



U 



W 



{n\y 



' -|^r^ 2n 

J2 (^k^...k2,A{Kki,---,Kk„)/\{Kk„+i,...,Kk2^) 

kl,...,k2n = l 

Xv^j.-, . . . v\i, Th-, . . . Tu , 
^ 2n 

J2 efc,...fc2„A(Kfc,, . . . , Kk„_,)A{Kk„, . 



(n-l)\(n + l)\ 



ki,...,k2„=l 



K, 



k2n 



xXk^ . . .Xk^^^Tk^ 



rk„ 



;i.lla) 



;i.iib) 



where 



A{K, 



.Kn) 



1 1 ■■■ 1 



;i.i2) 



Ki"-^ Kr^ ■ ■ ■ K~^ 

denotes the determinant of a Vandermonde matrix. On the symmetry axis, these formulas 
reduce to 



U 



W 



(n\ 



1 \n 2n 

J2 (^k^...k2u^{Kki,...,Kk^)A{Kk„+i,...,Kk2„) 

ki,...,k2„=l 

xXk,--- Xk^iKk, - z)--- {Kk„ - z) , 

■1 2n 

Y. ^k^...k2„MKk^, • • • , Kk„_,)A{Kk„, . 



in- l)!fn + l)! 



K, 



k2n 



kl,...,k2n=l 

xXk,--- Xk„_,{Kk, -z)--- {Kk„_, - z) 



;i.l3a) 



;i.i3b) 



which are consistent with the axis values of U and W being given by expressions ( |1.9a| ) and 
( |1.9b|) , respectively. 

II. THE iV = 2 SOLUTION 

The case n = 1 is well known to correspond to the Kerr-NUT spacetime. Therefore, we 
shall concentrate upon the next simplest case, n = 2, which was first considered by Kramer 
and Neugebauer,llil and where 

U = ~{K2 - K^){K, - K;){X,X2rir2 + X^X^r^u) 

+ {Ks - K,){K, - K2){X,X,nrs + X^X.r^n) 

-{K, - K^){K; - K2){X,X4rir, + X^X.r^r,) , (2.1a) 

W = -A{K2, Ks, K,)X,n + A(i^3, K,, K,)X2r2 

-A{K,, K,, K2)X,r, + A{K^, K^, K,)X,u ■ (2.1b) 

A simple calculation yields the following expressions for those complex constants that appear 
in the axis data: 

Uo= - {K, - K2){K^ - K^){X,X2 + X3X4) 
+ (i^i - K^){K2 - K,){X,Xs + X2X4) 

- {K, - K,){K2 - K,){X,X, + X2X3) , (2.2a) 
U^ = {K^ - K2){Ks - K,)[{K^ + 1^2)^1X2 + (K3 + K^)X^X,] 

- [K, - i^3)(i^2 - K^)[{K^ + K^)X,X^ + {K2 + ^4)^2X4] 

+ (i^i - K^){K2 - K,)[{K^ + K4)XiX4 + {K2 + K,)X2X,] , (2.2b) 

U2= - {K, - K2){Ks - K^){K,K2X,X2 + K^K^X^X^) 
+ (i^i - K^){K2 - K,){KiKsXiXs + 7^2^^4X2X4) 

- {K, - K^){K2 - Ks){K,K,XiX^ + ^^2^^3X2X3) , (2.2c) 
W, = A{K2, Ks, K^)Xi - A{K,, K^, K^)X2 
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+ A{K^, K2, K^)Xs - A{Ki, K2, Ks)X, , (2.2d) 

W2= - A{K2, i^s, K4)K,Xi + A{Ku K^, K,)K2X2 

- A{Ki, K2, K^)KsX^ + A{K,, K2, Ks)Ka4 ■ (2.2e) 

We were surprised how easy it was to solve these equations for the four complex constants 
Xa (a = 1, 2, 3, 4) in terms of the K's and the axis data. 

A. Determination of Xa {a = 1, 2, 3, 4) 

We begin by using the above expressions for Wi and W2 to express X4 in terms of the 
other X's. Then, we use the expressions for W2 and Ui to express X2 (and X4) in terms 
of Xi and X3. Next, we use the expressions for 14^2 and Uq to express X3 (as well as X2 
and X4) in terms of Xi. Finally, we use the expressions for W2 and U2 to solve for Xi (as 
well as X2, X3 and X4). The first and last equations are linear, while the second and third 
are, perhaps surprisingly, only quadratic. Choosing the roots of the quadratic equations 
judiciously, we obtain the following expressions for the X's: 

VX^ = [U2W2 + f/i(f/2W^i - U^W2)] + {K2 + K^ + K4){U2W^ - U1W2) 

- {K2K3 + K2K4 + K^K4)W2 - {K2K^K4)Wi , (2.3a) 

VX2 = [U2W2 + Ui{U2Wi - U1W2)] + (i^3 + K4 + Ki){U2Wi - U1W2) 

- {KsK^ + K^K^ + K4Kr)W2 - {K^K4K^)W^ , (2.3b) 
VX^ = [U2W2 + Ui{U2Wi - U1W2)] + {K4 + Ki + K2){U2Wi - U1W2) 

- {K4K1 + K4K2 + KiK2)W2 - {K4KiK2)Wi , (2.3c) 
VX4 = [U2W2 + Ui{U2Wi - U1W2)] + {Ki + K2 + K^){U2Wi - U1W2) 

- {K,K2 + K^K^ + K2K^)W2 - {K,K2K^)W, , (2.3d) 



where 



V := Wi{U2Wi - U1W2) + W^ . (2.4) 



In these expressions, the complex parameters Ua, Va, Wa {a = 1,2,3,4) have been rescaled 
so that Uq = 1. For each choice of the parameters Ka {a = 1,2, 3, 4), these equations assign 
values to the four complex parameters Xa (a = 1, 2, 3, 4). 

B. Determination of Ka (a = 1, 2, 3, 4) 

The requirement (|1.8| ) leads to the following four conditions (where Uq = 1): 

Ki + K2 + K^ + Ki = -2Re Ui , (2.5a) 

K^K2 + K^K^ + K^K^ + K^K^ + K^K^ + K^K^ = \Ui\^ - \Wi\^ + 2Re U2 , (2.5b) 

K2KsKi + K^KsK^ + K^K^K^ + K^K^K^^ = -2Re {U2U* - W2W*) , (2.5c) 

K1K2K3K4 = |[/2p - 11^2!' . (2.5d) 

Incidentally, these four relations are equivalent to the single relation 

\U{z, 0)p - \W{z, 0)|2 = \Uo\\K, - z){K2 - z){K, - z){K, - z) , (2.6) 

which one should also be able to deduce directly from the determinantal expressions for U 
and W, and which can be generalized in a natural way for n > 2.E3 

From the above expressions it is clear that each of the K^s satisfies the quartic equation 

= K^ + 2Re f/ifCf + (|f/i|2 - \Wi\^ + 2Re U2)K^ 

+2Re {U2U; - W2W*)Ka + (If/sP - 11^21') • (2.7) 

Moreover, because the coefficients of this quartic equation are real, the solutions Ka {a = 
1, 2, 3, 4) are real, or occur in complex conjugate pairs. 

Of course, using a translation along the z axis, we can always achieve Re t/i =0 or 
Ki + K2 + K3 + K4 = 0. Cardano's method of solving the quartic equation 

K^ - AKl - BKa + C = (2.8) 

then yields 



where 



Ki = -ik + ^2A - F + 2B/k\ , (2.9a) 

K2 = Uk- pA - P + 2B/k\ , (2.9b) 

Kz = \ l-k + ^2A - P - 2B/k\ , (2.9c) 

K4 = -l-k- ^2A - P - 25/A;| , (2.9d) 



^VT_ + -A , (2.10) 



T± := -^^(^' - 36C) + ^B' ± ^y^ , (2.11) 



and 



To := 815^ - 12B'^A{A^ - 3QC) - A8C{A^ - ACf . (2.12) 

In the vacuum case we have the identifications 

A = \Wi\^ - \Ui\^ - 2Re U2 , (2.13a) 

B = -2Re {U2U* - W2W^) , (2.13b) 

C=\U2\^-\W2\\ (2.13c) 

which can be generahzed easily for electrovac fields. 

An alternative way of expressing this result that facilitates a comparison with Manko et 
al. is 

K, = ^{k + {K+ + K.)} , (2.14a) 

K2 = ^{k-iK+ + K^)} , (2.14b) 

Ks = ^{-k + {K+-K^)} , (2.14c) 

K, = ^{-k-{K+-K^)} , (2.14d) 

where 



K±:=^^{A-k'^/2)±2d, (2.15) 

and 



d:=-^{A-k^/2f-{B/ky. (2.16) 

When the axis data happen to satisfy the relation 

B := -2Re {U2UI - W2WI) = , (2.17) 

one also has k = 0. Therefore, one must carefully evaluate the limit of our expressions for 
the K^s as -B — i> 0, noting especially that hraB^o B/k is finite, and lime^o'^ = VC- The 
reader can show that the result obtained this way is consistent with the fact that when 
5 = 0, the square of each Ka satisfies a quadratic equation, 

= iKl)'-AKl + C, (2.18) 

which can be solved directly. Without further loss of generality we can express the solution 
in the form 

Ki = -7^2 = ^(/€+ + /€_), (2.19a) 

K3 = -K4 = i(K+ - «:_) , (2.19b) 

where 



K^ ■= J\Wi\^ - |f/i|2 + 2{±d - Re f/2) (2.20) 



and 



d:=J\U2\^-\W2\^. (2.21) 



The reader can easily check that all equations ( |2.5a| ) through ( |2.5d| ) are satisfied by this 
solution. 



In Eqs. ( ^.9a| ) through ( [^.9cl| ) or Eqs. ( P^.14a| ) through ( |2.14dD we have expressed Ka 



[a 



1,2,3,4) explicitly in terms of the axis data. Either of these sets of expressions can be 
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substituted into Eqs. ( p.3a|) through ( |2.3d| ) to obtain Xa (a = 1, 2, 3, 4) exphcitly in terms of 
the axis data. The complex potential £ is then given by Eq. ( |1.3| ) with U and W rendered 
by Eqs. ( p. la} ) and ( |2.1b| ), respectively. 



III. THE SPACETIME METRIC 

In principle the metric, which is usually written in the form 

ds^ = f-^[e^\dz^ + dp^) + p^d(j)^] - f{dt - ud(Pf , (3.1) 

can be constructed once the complex potential S is known. The field / := —gu = Re £^ is 
given by 

Thus, the infinite redshift surface corresponds to 

|[/|2 = \W\'^ {Infinite Redshift Surface}, (3.3) 

while it can be shown that there is a curvature singularity whenever 

U + W = {Curvature Singularity}. (3.4) 

The fields u and 7 can be determined, up to respective integration constants, by solving the 
differential equations 

i^P = -pf'^Xz , (3.5a) 

<^z = pf'^Xp , (3.5b) 

^, = p{s,s;-SiS;}/{2f)\ (3.5c) 

^. = p{s,s; + SiS;}/{2f)\ (3.5d) 

where x •= I^^ ^■ 

While, in principle, this is straightforward, in practice it is extraordinarily tedious to 
calculate the field u in this way. Even for n = 2 the number of terms one encounters is 
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enormous. Modern solution-generating techniques usually provide some alternative method 
for determining oo. We are most familiar with our own homogeneous Hilbert problem (HHP) 
formulation, in which the iJ-potential of Kinnersley and Chitre plays a key role. This is a 
2x2 matrix potential, the negative of the real part of which is just the metric block 



h 



' f-'p' - f^' fco ^ 



V 



(3.6) 



/ 



while the lower right element of the matrix H is the complex £^-potential. Now, the HHP 
yields not just a formula for S, but rather a formula for H; namely. 



H = H^^^ + x+{o)n 



(3.7) 



where, for Minkowski space. 



/ .2 n\ 



H(o) 



and, generally. 



n-.-- 



-p' 

-2iz 1 



fo A 



v-'»/ 



(3i 



(3.9) 



The matrix X^ (t) is one of the 2x2 matrices involved directly in the HHP. The S potential 
is given by 



s = ncn = s^'^ + )Kif: 



cc 



(3.10) 



while 



fuj = -Re H^^ = /(o)cj(°) - Im X+(0)^^ , 



(3.11) 



where LL and LR refer to lower left and lower right matrix elements, respectively. The point 
is that one can evaluate X^{0)^^ almost as easily as X^{0)^^. The field u is determined 
much more easily this way than by integrating the differential equations for u. 
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For the general n = 2 solution of the Neugebauer family, the procedure we have outlined 
yields a result of the form 

f{u;-u;o) = -lm [jf^) , (3.12) 

where 

Af = {K, - K2){Ks - K,) 

x{K^ + K2-K,- K,){Q^Q2 - QsQ,) 
-{Ki - Ks){K2 - K,) 

x{Ki + Ks-K2-K,){QiQ,-Q2Q4) (3.13) 

+ {K, - K,){K2 - Ks) 

x{Ki + K,-K2- Ks){QiQ4 - Q2Q3) , 
Qk:=t[Xkrk + iKk-z)], (3.14) 

and Uq is a real constant, the value of which is determined by the HHP in such a way that 
cj = on the axis. 

Finally, one finds that, for all n = 2 solutions, the field 7 is given by the simple expression 

e^T = (|t/|2 - \W\^)/{\Uo\W2r3n) , (3.15) 

where the constant has been chosen so that 7 = on the symmetry axis.ll^ 
In particular, if we specialize to the case 

Ui = -ia , U2 = b , Wi = m , W2 = , {a,b,m real) (3.16) 

where a^ < m^, we find that 

X, = X, = ~^±±^ ,X2 = X, = ^±^^ . (3.17) 

m m 

Substituting these values into our general expressions for U, W, and J\f, and using the 
notation of Manko et al., 
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i?+ := r2 , R^ := ri , r+ := r^ , r_ := rs 



(3.18) 



we obtain 



[/ = % ((m^ - 2a2)(i?_r„ + R+r+) + 2iaK+(R^r_ - R+rM 

+/t^(i?_r+ + R+r_) - 4b{R+R_ + r+r_) , 
W = ^^^^^ {(m^ - a'^){r+ + r.- R+- i?_) + K+/t_(i?+ + i?_ + r+ + r„ 
+m[(fi;+ + K-){r- - r+) + (k+ - K-){R+ - R-)]} , 



m 



(3.19a) 



(3.19b) 



the vacuum specialization of the most recent solution published by Manko et al. Our way 
of expressing the field cj is a good deal simpler than is theirs; namely, 



f{uj - Uq) = -Im 



U + W 



where loq = —2a, and 



AT = €/€+/€_ {k4Q^Q^ - Q^Q^) - k+{Q^Qa - QeQ3)} 



where 



Qi =i 
Q2 = i 
Qs = i 
Qi = i 



K+ +ia\ „ 1/ , ' 

m / 2 

Kj^^ — ia\ ^ 1 . , 

m / 2 

K+ + ia\ 1 

m / 2 

K^ — ia\ 1 

m / 2 



(3.20) 



(3.21) 



(3.22a) 
(3.22b) 
(3.22c) 
(3.22d) 



IV. FUTURE EXTENSIONS 

In this paper we have succeeded in expressing all the n = 2 members of the Neugebauer 
family of solutions of the vacuum field equations in terms of data prescribed on the sym- 
metry axis, which in turn are easily related to the multipole moments of the source of the 
gravitational field.B'ty One first evaluates the K^s using Eqs. ( p.9a|) -( p.9d| ) and the X's using 
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Eqs. (|2.3a|) - (|2.3d|) . Then one evaluates U and W using Eqs. (p.la|) and (|2.1bD , respectively, 
and fuj using Eqs. ( |3.12| ), ( p.l3| ) and ( |3.14| ). Of course, 7 is given by Eq. ( p.l5[ ). 



In a future paper we shall include electromagnetic fields, where the complex potentials 
have the forms 

U-W V 

^ - ^TTvv ' ® " ll+w ' ^^-^^ 

and U is homogeneous of degree n while V and W are homogeneous of degree n — 1 in 
ri, . . . , r2n- Electrovac fields of the type in which we shall be interested have been generated 
by at least two techniques, one due to Alekseevllil and the other due to Cosgrove.ll3 We are 
more familiar with Cosgrove's approach, which, in its usual formulation, produces directly 
only the charged Kerr metric with a^ + e^ > m'^. The complexified Cosgrove transforma- 
tion, in which the group SU{2,1) is replaced by SL{3,C) produces a family of complex 
spacetimes,El the real cross sections of which are the electrovac spacetimes we shall study. 

After considering the electrovac extension, we shall pass on to the case n > 2. For all 
values of n, the K^s will satisfy the equation 

2n 
a=l 

on the symmetry axis.ll3 Except in special cases, it will not be possible to express the K^s 
as algebraic expressions in the axis data Ua, K, Wa {a = 1, . . . ,n), because, in general, the 
K^s will be solutions of an algebraic equation of degree 2n. On the other hand, it may 
be possible to express the complex potentials S and $ explicitly in terms of the axis data 
and the K's, with the latter parameters determined numerically from the axis data or the 
multipole moments. 

Long ago Neugebauer gave the 2n x 2n determinants for U and W, but we are not aware 
of any similar formula for feu. Likely general formsli^ for J\f and 7 can be guessed from the 
expressions we have displayed for the case n = 2: 

1 2n 

^=r-n^ Y. ^k,...k2nMKk„...,Kk„)A{Kk„+^,...,Kk^J 

X [{Kk, + ... + KkJ- {Kk^^, + ... + Kk,J]Qk, ...Qk„, (4.3a) 
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and 

2n 

exp(27) = (|f/p - |l^n/(|f/oP n ^^) • (4-3b) 

fc=i 

We regard it as extremely unlikely that any amount of study of the particular instances 

considered by Manko et al. would have allowed one to guess a plausible generally applicable 

form for uo. That is the principal advantage to considering the general n = 2 case rather 

than particular examples, no matter how interesting those particular examples might be 

with respect to potential physical applications. 

It should also be mentioned that, following Neugebauer, one can also generalize our work 

by using an arbitrary Weyl metric as the seed metric instead of Minkowski space. One can 

conceive of potential physical applications for some of these metrics too. 
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